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Understanding interactions between membranes requires measurements on well-controlled systems 
' close to natural conditions, in which fluctuations play an important role. We have determined, by 

grazing incidence X-ray scattering, the interaction potential between two lipid bilayers, one adsorbed 
■ on a solid surface and the other floating close by. We find that interactions in this highly hydrated 

' model system are two orders of magnitude softer than in previously reported work on multilayer 

stacks. This is attributed to the weak electrostatic repulsion due to the small fraction of ionized 
^\ . lipids in supported bilayers with a lower number of defects. Our data are consistent with the 

^SJ ' PoissonBoltzmann theory, in the regime where repulsion is dominated by the entropy of counter 

ions. We also have unique access to very weak entropic repulsion potentials, which aUowed us to 
discriminate between the various models proposed in the literature. We further demonstrate that the 
' interaction potential between supported bilayers can be tuned at will by applying osmotic pressure, 

O ' providing a way to manipulate these model membranes, thus considerably enlarging the range of 

biological or physical problems that can be addressed. 
-4— • , 
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O ■ I- INTRODUCTION 

Supported lipid bilayers offer a unique configuration whereby a single bilayer, accessible to other molecules such 
as, for example, proteins, peptides, or DNA, is supported on a solid substrate. Beyond their interest for biosensor 
^ ' technology, the access they give to a flat immobilized membrane makes them highly relevant for fundamental studies 
in biophysics and membrane biology [l|, [^j. particular, they provide a unique way to finely characterize the inter- 
actions between membranes and their environment, which are not only crucial for membrane fusion and trafficking, 
endocytosis, and exocytosis 0, 3, but also fascinating from the physical point of view. 

Membranes indeed exhibit extremely complex interactions with their environment, in which both molecular-scale 
enthalpic and fluctuation-related entropic contributions are inextricably involved. In particular, the effect of confine- 
, ment has been now discussed for 40 years without a definitive answer being found. Helfrich first realized that, in 
addition to the direct electrostatic, van der Waals, and hydration forces Q,the long-range effective steric interaction 
generated by the thermal fluctuations of confined flexible membranes is an essential contribution to the total free 
energy of interaction [sj. Pure hard wall interaction (hard confinement) was first considered in [11,01 but is not a 
realistic description of real systems, and especially not of living ones. Confinement by a "soft" potential was treated 
I . either by using self-consistent methods leading to effective exponentially decaying potentials Q-Q, or by estimating 
average values within a full statistical mechanics approach [l3]. Which functional form should be used to describe 
entropic repulsion in real experimental situations, however, remains an open question. 

Although the surface force apparatus can be used to precisely determine the direct part of the potential (hydration 



and van der Waals contributions) [ll|,the entropic repulsion can only be studied by using scattering techniques. 
Combination with osmotic pressure measurements allowed in particular the determination of the compressibility 
B = d^J-/dd^, where F is the system free energy and dyj the interlayer water thickness [ll|, [13 ■ Although some 
agreement was found with the soft potential of Ref.Q, the experimental decay lengths found in [l^ were greater 
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FIG. SI: Schematic view of the experimental setup for specular and ofF-specular reflectivity. The grazing and scattered 
wavevectors (resp. angles of incidence) are ki„ and k^c (resp. 6in and 9sc)- q is the wavevector transfer. 



than twice the value predicted by theory, ^ 0.2 nm. Moreover, inconsistencies between the temperature dependence 
of B and the observation of an unbinding transition pointed to the role of static defects in multilayers, which would 
dramatically affect bilayer interactions |14| . In the work described here, we determine the interaction potential between 
supported bilayers. These bilayers can be almost defectless but contain much less material than multilayers and could 
not be studied using diffuse scattering until recently [15 1. 



II. RESULTS AND DISCUSSION 



Two kinds of supported bilayers were investigated in this study, both consisting of two bilayers (see Fig. |T]). The 
first type, called "double bilayers" in the following, consists of two bilayers of DSPC, whereas in the other, called 
"OTS-bilayer" , the first monolayer close to the substrate is replaced by an octadccyl-trichlorosilanc grafted layer 
(see Materials and Methods for details). In both cases, the second bilayer is free to fluctuate in the potential of the 
first bilayer and of the substrate. A combined fit of specular and off-specular data is performed in order to increase 
sensitivity, using a model taking into account the static and thermal roughness from both bilayers (see Ref. [l6j 
and Materials and Methods). From the fits presented in Fig. [S2l we obtain structural parameters, in particular the 
bilayer-bilayer distance and the interlayer water thickness, but also the interaction potential second derivative, and 
the bilayer tensions and bending rigidities (Fig. [S3]). Electron density profiles for OTS or double bilayers can be 
found in Ref. (l6j (an example is given in Supporting Information). 

In all cases, the best fit values for the fioating membrane tension is 0.3 ± 0.2 mN/m. The bending rigidity decreases 
from (250 ± 50)fcsT in the gel phase (T = 42.9°C) to (50 ± 20)kBT in the fluid phase (T = 52.9°C) in agreement 
with previously reported values [isj . The static, substrate-induced roughness of both membranes is always less than 
0.3 nm and remains constant, as well as the thermal roughness ath,i of the flrst, adsorbed membrane, which is on the 
order of 0.4 to 0.5 nm (Fig. [S4] inset). The thermal roughness ath,2 of the second bilayer is larger than that of the 
first bilayer, justifying the denomination "fioating bilayer", in good agreement with previous experiments [l6l.[l7j. 
A zoom of the off-specular reflectivity in the region where it is most sensitive to the interaction potential is shown in 
Fig. [S2]for different temperatures. It is important to note here that the second derivative of the interbilaycr interaction 
potential is directly linked to the depth of the minimum in the diffuse scattering curve around qz ~ 1.0 — 1.5 nm~^, 
without much coupling to the other parameters. Similarly, the interlayer water thickness is strongly correlated to the 
Qz position of that minimum. Hence, it can be seen directly in the inset of Fig. IS2I that the interaction potential 
becomes weaker (the minimum is less pronounced) when the interlayer water thickness increases (left shift of the 
minimum) at higher temperatures. 

The second derivative of the interbilaycr potential obtained by fitting the experimental data is represented 
as a function of the interlayer water thickness in Fig. IS3[ where our data are compared with values obtained 
by Pctrachc ct al. for Egg PC multilayers [l^. Remarkably, our samples arc more hydrated than multilayers {dw 
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FIG. S2: OfT-specular reflectivity from silicon substrate (o) and an OTS bilayer at T=42.9°C (□) as a function of q^. Continuous 
lines represent best fits. Top inset: off-specular reflectivity as a function of qz zoomed in the region where it is most sensitive 
to the potential at T=42.9°C (□), 49.7°C (■), 51.5°C (o), 52.0°C (•), 52.9°C (A). Note the shift in minimum and decrease in 
contrast with increasing temperature. Bottom inset: specular reflectivity (continuous line), "true" specular reflectivity (dashed 
line), and diffuse scattering in the specular direction = (dotted line) calculated using the model of Ref. [l^ using the flt 
parameters and experimental resolution. 



is 0.1 to 0.5 nm larger) and interact via a softer interaction potential (//"g is smaller). Note that in this analysis, 
the interaction potential is not only strongly constrained via its second derivative U12, but also via dw which fixes 
the position of its minimum, is in turn reported as a function of the second bilayer thermal roughness ath.2 in 
Fig. IS4I which shows a strong correlation between the two parameters, demonstrating that there is a large cntropic 
contribution to the repulsion as expected. 

The accuracy of our data, in particular for large separations, allows for a precise assessment of the repulsive part of 
the potential. In particular, it allows for a test of the different functional forms which have been used in the literature 
to model the entropic part of the potential. With this aim, we first accurately calculated the attractive part of 
the potential using the Lifshitz theory [l^ , carefully modelling the silicon-silicon oxide- watcr-lipid bilayer- water- lipid 
bilayer stack (for more details see Supporting Information). For du, < 3 nm, a good approximation to the van der 
Waals interaction energy is ?7vdw = ~H/12Tr{dw + Sdhcad)^, with H = 5.3 x 10~^^ J and dhoad the headgroup layer 
thickness (0.4 — 0.8 nm) in good agreement with previous work [l9l. [20|. 

We further modeled hydration forces using a classical exponential decay, Uhyd = PhZh exp {—z/zh), with = 1 — 5 x 
10'' Pa the hydration pressure and Zh = 0.16 — 0.2 nm the hydration length p^ . 

The renormalization of the microscopic interaction potential by the thermal fluctuations is a complex problem of 



modern statistical physics. According to Helfrich [21[, the membrane free energy is the sum of the microscopic 
potential and of the entropy cost of confining the flexible membrane, which results in an effective potential controlling 
bilayer position. The effective potential, average bilayer position, and fluctuation amplitude are thus coupled quantities 
which must be sclf-consistently determined. 

As mentioned in the introduction, Podgornik and Parsegian have extended Helfrich's approach [2l| to take into 
account hydration repulsion and van der Waals attraction in the so-called "soft" potential Q , and the self-consistent 
approaches of Rcfs. [13, 1111 allow the calculation of average position, rms roughness and mean effective potential 
curvature. All these theories are based on a quadratic approximation of the interaction potential, either symmetric 
[2H or not [13, ^'^y case, a non-symmetric case, like ours, can always be mapped to a symmetric case by 

identifying the strength of the quadratic potential with the second derivative of the asymmetric potential. As an 
example, the Helfrich potential has been shown to correctly describe a supported bilayer interacting with a single 
hard wall in the limit of small fluctuations (see [loj and Supporting Information) . 
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FIG. S3: Interaction potential second derivative U12 as a function of the interlayer water thickness d^. (□): data from Petrache 
et al for Egg PC in fluid phase. All other data from this work using DSPC: fluid phase (•); fluid phase prior to (•) and 
after(o) the addition of salt (T=58°C; £d ~ 0.5 nm and 0.4 nm); gel phase prior to (a) and after (A) the addition of salt 
(T=40°C; Id = 0.4 nm); (♦) Gel Phase with applied osmotic pressure (II = 10 ±4 and 450 ±250 kPa). Red solid (respectively 
dashed) lines: soft-confinement potential plus electrostatic contribution in the Ideal Gas limit (IG) for Fluid (respectively 
Gel) phase. Blue solid (respectively dashed) lines: soft-confinement potential plus electrostatic contribution in the Debye- 
Hiickel limit (DH) {£d = 0.3 nm) for Fluid (respectively Gel) phase. Green dashed-dotted line: self-consistent model in 
fluid phase (k = 50fcsT). The ellipses show the region of equilibrium states without applied pressure for Helfrich, hydration 
and van der Waals forces (grey area); soft-confinement, hydration and van der Waals forces (dashed area); soft-confinement, 
hydration and van der Waals forces plus electrostatic interaction in the ideal gas limit (empty ellipse). The difi'erent ellipses 
were obtained by varying Ph, Zh, dhcad and k within the limits indicated in text. Inset: osmotic pressure 11 (interaction potential 
first derivative U12) as a function of the interlayer water thickness dw (□): Open symbols, data from Petrache et al [T^] for 
Egg PC (□), DMPC (0) and DPPC (A) in fluid phase. (♦) data from this work using DSPC in gel phase with applied osmotic 
pressure. Red dashed line: soft-confinement potential [Q] with electrostatic contribution in the Ideal Gas limit for the Gel 
phase. Blue dashed line: soft-confinement potential plus electrostatic contribution in the Debye-Hiickel limit (DH) (^d=0.3 
nm) for the Gel phase. Green dashed-dotted line: self-consistent model [T^] in gel phase (k = 250A;bT). Black dotted line: 
microscopic potential plus electrostatic contribution without any entropic contribution. The same set of parameters has been 
used in main figure and in the inset. 



We first considered hard wall repulsion, with Uuei = ch / ^{kBT / z)^ [2l| per unit area and Ch ranging from 0.08 to 
0.2 [lol.[23|. The corresponding d^j equilibrium values at zero applied pressure are given by the dark grey area on Fig. 
[S3] and show poor agreement with the data. As expected, the hard- wall potential also fails to describe the vs (Jth,2 
curves (Fig. [S4|) . Accordingly, the classical "Helfrich Ansatz" d^, oc ath does not apply (Fig. [S4)) . We also note that 
simply shifting the zero of the potential to account for softness does not help and that the potential of Ref. Q taking 
membrane tension into account cannot be distinguished from the hard-wall potential for realistic tension values. The 
"soft" potential of Podgornik and Parsegian Q J7soft = nksT / 16 y^I\jKjzhexp{—z/nzh) leads to a slightly better 
agreement but still predicts U12 values one to two orders of magnitude larger than those observed experimentally 
(Fig. IS3I) a nd also underestimates the values of d^ and ath.2 data (Fig. IS4|) . Finally, the self consistent theory 
of Ref. [lO[, which in principle allows one to calculate the mean d^, [/"j or (Tth,2 more satisfactorily than effective 
potential theories, is indeed in good agreement with the d^^ vs ath.2 data, but it strongly overestimates c?^, ath and 
U12, probably because of a bad sampling of the most confined microstates. 

As our bilayers weakly interact with the substrate and can be very close to detachment, which was sometimes 
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observed for temperatures > 60° C, a very small but long-range repulsive contribution to the potential would in fact 
be enough to shift the equilibrium position to higher and lower [/{'j. Such electrostatic interactions are investigated 
in Ref.fiS]; but always ignored in scattering studies on zwitterionic lipids [l^l- Phosphatidylcholines present indeed 
pKa values of 2.7 and 11 [1^ and bear a positive charge density a ~ 0.001 e~/nm^ at the experimental pH = 5.5. 
This small amount of charges, necessarily present due to the amphoteric character of the phosphatidylcholine group, 
leads to a weak electrostatic repulsion which was recently shown sufficient to prevent vesicles from adhering [26| . 




FIG. S4: dm as a function of ath.2- Symbols are the same as for figure [S3l Dashed-dotted line (blue online): "Helfrich Ansatz" 
dw = ^/GcFth 2- Dotted lines (blue online): Helfrich confinement, with a 0.3 nm shift for the lower curve as explained in text. 
Dashed line (red online): Soft potential (Ref. [^) without electrostatic contribution. Solid line (red online): Soft potential 
(Ref. [3|) with electrostatic contribution. Green dashed-dotted line: self-consistent model [lo| with electrostatic contribution. 
Inset: static roughness aat,i (open symbols) and thermal roughness ath,i (closed symbols) of the first strongly adsorbed bilayer 
(i=l, square, red online) and of the floating bilayer (i=2, triangle, blue online). 



In our experiments, surface charge densities are small and the mean-field Poisson-Boltzmann (PB) theory is 
expected to appropriately describe the system. For DSPC supported bilayers in ultra-pure water, the Gouy-Chapman 
length related to the surface charge density cr is 10^ — 10^ nm, the Debye-Hiickel length which describes 

screening is in ~ 200 — 500 nm and ^2 — 3 nm. In this regime, the Ideal Gas (IG) limit of the mean-field 
Poisson-Boltzmann theory, where the interaction potential reads Uei ~ 2kBT(j / e hi{d.uX, should apply rather than the 
Debye-Hiickel limit (DH) [23| ■ Adding this contribution to the soft potential of Ref. [9[ gives a perfect description of 
our data for both representations vs dw (Fig. [S3)) and d^ vs ath,2 (Fig. [S4)) . Adding the electrostatic contribution 
extends the curves previously obtained to larger <Tth.2 and d^ without significantly modifying their low d^, part. 
Therefore, it cannot significantly improve the agreement for the other forms of the potential (Fig. IS3[) . Accordingly, 
the charge density is very robust and docs not depend on the other parameters. The validity of the soft effective 
potential in the presence of electrostatic interactions could nevertheless be questioned as the soft potential was 
explicitely constructed for hydration and van der Waals forces only. Though the infiuence of electrostatic interaction 
on entropically induced repulsive forces has not been investigated in detail, these are very long range compared to 
van der Waals and hydration interactions, and should be only marginally renormalized. Moreover, it was shown in 
Ref. [l^l, that in the absence of van der Waals and hydration forces, the renormalization of electrostatic interactions 
is negligible in the limit where the in-plane electrostatic correlation length ^ (50-100 nm here) is larger than d^,. 
In order to further check the effect of electrostatic interactions, NaCl was added to a double bilayer. This is expected 
to increase screening and is easier to analyse than changing the pH to change the lipids degree of ionization. We 
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prepared two solutions with in = 0.45 nm and Id = 0.3 nm (see Materials and Methods). The results are reported 
in Fig. [53] and [54] A strong decrease in the interlayer water thickness and in the thermal roughness ath,2 is 
observed, as well as a large increase in the interbilayer potential C/{2, in good agreement with the strong screening of 
the electrostatic potential. 

A further proof of the good control we have over the interactions between supported bilayers is provided in the inset 
of Fig. IS3I where the effect of osmotic pressure is shown. Osmotic pressure was applied using PVP (see Materials and 
Methods for details), and here again, the agreement with the theoretical model and with previous experiments on 
multilamelar systems is perfect. This is a first demonstration that the interbilayer potential of supported bilayers can 
be tuned using osmotic pressure, allowing us to extend the measurements towards smaller interlayer water thicknesses 
and to bridge the gap with multilayer studies. 

III. CONCLUSION 

In this article. X-ray off-specular scattering measurements of the interaction potential between two bilayers adsorbed 
on a solid substrate are shown to lead to results presenting unprecedented sensitivity, illustrated by the necessity of 
taking into account the very weak electrostatic repulsion between almost neutral bilayers and the possibility of 
discriminating between different entropic and electrostatic potentials. These results show that supported bilayers 
are significantly more hydrated and therefore exhibit more intrinsic properties than the usually studied multilayers, 
possibly owing to defects in the latter. This opens up a wide range of possibilities for understanding unbinding or 
investigating the effect of various biological molecules on interaction and adhesion between membranes. 

IV. MATERIALS AND METHODS 

The supported bilayers were prepared by depositing two bilayers on ultra-flat silicon substrates (SESO, France), 
where the first, more strongly adsorbed, bilayer is either a bilayer of L~a I,2-distearoyl-sn-glycero-3-phosphocholine 
(DPSC, Avanti Polar Lipids, Lancaster, Alabama) made by a combination of classical Langmuir-Blodgctt (LB) and 
Langmuir-Schaefer (LS) depositions (vertical sample) [l^ or a mixed octadec ytri chlorosilane (OTS) - lipid bilayer 
(OTS bilayers), where the OTS layer is chemically grafted on the substrate [2^. A second, "fioating" bilayer is 
then prepared by a LB deposition, followed by a LS deposition. The first bilayer serves both as a spacer to reduce 
the interaction between the floating bilayer and the substrate and keep it free to fluctuate, and to investigate 
bilayer- bilayer interactions. The sample is then inserted into a PTFE sample cell with SO/^m thick windows. 
The sample cell is tightly closed and transferred to an alumina box, thcrmalized by a water circulation, first at 
25°C, then heated by steps, with a feedback on the temperature measured inside the sample cell by a PTIOO resistance. 

Experiments were performed using ultra pure water (18.2 Mn.cm), obtained from a Millipore purification system. 
Dissolution of CO2 in water leads to Debye-Length values around ~200 nm, ten times smaller that the 960 nm expected 
for such samples [s^]. By adding Sodium Chloride c=0.5 and 1 M, we obtained solutions with Dcbye length equal to 
0.45 and 0.3 nm respectively. Osmotic pressure was applied using polyvinylpyrrolidone (PVP) of average molecular 
weight 40 000 (Sigma Aldrich, St. Louis, MO) mixed with milli-Q water. Solutions of 4 and 30 % of PVP (w/w) 
were prepared and homogenized with magnetic stirrings overnight. The bulk water was then carefully replaced by the 
PVP/water solutions with syringes, taking care not to expose the bilayers to air. The values of the osmotic pressures 
are deduced from the PVP concentrations as calculated by Vink [3l| and reported in [l^, leading to n= I4±2 kPa 
and 450±250 kPa. 

Specular and off-spccular reflectivities were recorded using the procedure of Ref. [isj. The experiments reported 
here used a 27 keV x-ray beam (wavelength a=0.0459 nm) at the CRG-IF beamline of the European Synchrotron 
Radiation Facility (ESRF). The scattering geometry is described in Fig. H] The monochromatic incident beam was 
first extracted from the polychromatic beam using a two-crystal Si(lll) monochromator. Higher harmonics were 
eliminated using a W coated glass mirror, also used for focusing. In all experiments, the incident beam was 500 /^m x 
18 /xm (H X V). The refiected intensity was defined using a 20 mmx 200 ^m (H x V) slit at 210 mm from the sample 
and a 20 mm x 200 /^m (H x V) slit at 815 mm from the sample and recorded using a Nal(Tl) scintillator. Specular 
reflectivity was obtained by rocking the sample for each angle of incidence {q^ scans for approximately constant ) in 
order to subtract the background. Off-specular reflectivity was measured at a constant grazing angle of incidence of 
0.7 mrad below the critical angle of total external reflection at the Si- water interface (0.83 mrad), leading to variation 
of both q^ and (see Fig. |T|). 
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Optimization of slit widths allowed us to extend the in-plane wavevector transfer range gy by one order of magnitude 
from 2x10'^ m^i to 2xio^ m^^compared to the experiments of Ref. [l^. Our experiment is thus sensitive to the off- 
specular diffusion by both bilayers and not only the more strongly fluctuating one, allowing for a precise determination 
of interaction potentials. 

The differential scattering cross-section (power scattered per unit solid angle per unit incident flux) can be written 



PSi 



(1) 



where re=2.8i8xio"^'' m, t(ei„) and tie^c) are the Fresnel transmission coefhcients between water and silicon, for the 
grazing angle of incidence ei„ and for the scattering angle e,^. The coefficient t(6ii„) represents a good approximation to 
the actual field at the interface while t{e^^) describes how the scattered field propagates to the detector. Eq. ([T|) has 
to be be multiplied by the incident flux and (numerically) integrated over the detector solid angle to get the scattered 
intensity. The three terms in between the square brackets describe the surface roughness and bilayers 1 and 2 static 
roughness and thermal fluctuations respectively. Spi{q^) (i=l,2) is the Fourier transform of the i-th bilayer (located 
at zi(r,|)) electron density profiles (form factors), which are described using the so-called IG-hybrid model ficl. [33j. 
Expanding the square modulus in Eq. ([1]) , we get self- and cross- height-height correlation functions of the substrate 
and bilayers, where the cross-correlations are sensitive to the interaction potentials. 

We describe the substrate correlation function using a self-affine correlation function [s^] ■ Static and thermal corre- 
lation functions used in Eq. ([T]) arc derived in detail in Ref. [l6j using the free energy: 



1 ^ 

oE |5^(9||)P + C/;^5^('?||)5s(-9||)] -C/f2^l(9||)^2(-<7||), 



(2) 



with a(g||)=c/,"-i-7igj|+Ki(jJ}, whcrc 7i and Kt are respectively the tension and the bending modulus of the i-th bilayer. 
u"^ and u['2 are second derivatives of the effective interaction potential between the substrate and a bilayer and 
between bilayers. The linear response theory of Swain and Andelman [ssj was extended to double bilayers in order 
to describe the static coupling of the bilayers to the substrate and the thermal correlation functions were derived by 
diagonalizing .f,, applying the equipartition of energy and Fourier transforming. 
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SUPPORTING INFORMATION 



I. FIT PARAMETERS 

As described in the manuscript, we perform combined fits of experimental specular and off-specular data. Best fit 
main parameters for elastic constants (72, K2), interaction potential second derivative {U'j^j_^ a/^)' static {a'st,2) Stiid 
thermic {(yth,2) roughnesses of the floating bilayer are given in Tables [ST] and [S2l A typical Electron Density Profile 
(EDP) for a double supported bilayer is also reported as Fig. [55l 

A. Influence of temperature and osmotic pressure 





Influence of temperature 


Osmotic pressure 


42.9°C 49.7°C 51.5°C 52.0°C 52.9°C 


PVP 4% PVP 30% 


D2,H20 [A] 


25.0±0.2 26.5±0.2 26.9±0.2 28.1±0.2 29.7±0.2 


18.0 ±1 14.0 ±1 


K2 [ksT] 
72 [mN/m] 

(72, St [A] 
cr2,th [A] 


280±50 300±50 300±50 60±20 50±20 
0.3±0.2 0.3±0.2 0.4±0.2 0.4±0.2 0.5±0.2 
2.0±0.5 2.5±0.5 3.5±0.5 3.0±0.5 3.0±0.5 
6.3±0.7 8.6±0.9 10.7±0.8 12.1±2 13.2±3 


200±50 500±50 
l.llzb 0.3 0.3±0.1 
3.5±0.5 2.5±0.5 
3.8±1 2.8±1 




12.05±0.3 11.8±0.3 11.6±0.3 lf.7±0.3 11.5±0.3 


13.85±0.1 14.5±0.1 



TABLE SI: Temperature evolution and osmotic pressure effect. 



B. Influence of salt (NaCl) 





Fluid phase 


Gel phase 


Before salt 


Zd=0.5 nm 


?D=0.4 nm 


Before salt 


Zi5=0.4 nm 


D2,H,0 [A] 


28.0 ±1 


24.6 ± 1 


20.3 ± 1 


24.3 ±1 


22.2 ± 1 


«:2 [ksT] 
72 [mN/m] 
0-2, St [A] 

cr2,th [A] 


80 ±50 
0.80 ± 0.1 
0.5±0.5 
7.5 ±1 


473 ±50 
0.15±0.1 
1.0 ±0.5 
4.9 ±1 


232 ±50 
0.79 ± 0.1 
2.5 ±0.5 
4.1 ± 1 


75 ±10 
0.69 ± 0.1 
1.0 ±0.5 
8.0 ±1 


175 ± 10 
0.51±0.1 
1.0 ± 0.5 
5.8 ±1 




12.1 ±0.1 


12.9 ±0.1 


13.3 ± 0.1 


12.2 ±0.3 


13.0 ±0.1 



TABLE S2: Effect of salt screening in fluid and gel phase. 
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C. Example of Electron Density Profile 



0.25 



0.00 



Before adding PVP 
After adding PVP solution (4 %) 




4 8 

z (nm) 



A 



12 



FIG. S5: Example of Electron Density Profiles showing the effect of osmotic pressure. The upper curve has been shifted of 0.4 
e~/A^ for clarity. 



II. INTERACTIONS : LIFSHITZ APPROACH 



The Lifshitz theory describes the electromagnetic interactions between two continuous media in term of fluctuations 
of the electromagnetic field, summing over all frequencies. This approach is more general than the simple sum of 
pairwise interactions between atoms. For a review see [l|. 

We assume a realistic model of supported double bilayers (see Fig. IS6[) . with two semi-infinite media (water and 
silicon), two layers representing the two lipid bilayers (thicknesses &i and 62) and two water layers (thicknesses li and 
Z2). Lipid heads are highly hydratcd and can be included in water layer. So in our case ^2 = c^iij + 2dhead where dhead 
is the lipid head thickness. Using these model, Lifshitz theory leads to [l|: 



knT 



87r/i — 



E 



xln [(1 - Asj,^i (/i,&i) A.u,2,^^ (62)6 

X {l-As^,n^, {ll,bi)A (62)e-")]da;, 



(3) 



where A and A depend on the different geometric parameters {bi, U) and on dielectric constants of each medium. 
The full expression can be found in reference [l|. The key point in such a calculation is to accurately describe the 
dielectric permittivity frequency dependence. Following [2[, (cj) is given by : 



£4(1.^) 



iu)/ujr ^ 1 - w^/w^ + i7pa;/w^ 



(4) 



The uJr,p are the absorption frequencies, Cr,p the absorption strength and the ^r,p are the damping factors. The first 
sum corresponds to microwave relaxation and the second one to experimental absorption peaks, usually corresponding 
to Infrared and Ultraviolet frequencies (see Table IS3P . As shown by Ninham and Parsegian Q , in biologically relevant 
cases like lipid-water systems, it is of high importance to use the total dielectric data from component substances 
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FIG. S6: Schematic representation of a supported double bilayer for the van der Waals potential calculation. 





Microwave 


Infrared 


Ultraviolet 






IjJMW 


Cmw 




CiR 


'yiR 




Cuv 


7c/ V 






[10" rad.s"'] 




[lO^rad.s-^] 




[10'^ rad.s"'] 


[10^*^ rad.s-^] 




[10^^ rad.s"'] 


Silicon 


11.6 






0.345 


0.043 





0.503 


10.448 













0.535 


0.050 



















0.884 


0.059 











Water 


80.1 


1.083 


75.3 


0.314 


1.4635 


2.29 


1.2593 


0.0392 


0.774 










1.047 


0.7368 


5.78 


1.5172 


0.0570 


1.335 










1.396 


0.1526 


4.22 


1.7296 


0.0923 


2.336 










3.065 


0.0136 


3.81 


1.9724 


0.1556 


3.110 










6.450 


0.0751 


8.54 


2.2606 


0.1522 


4.491 
















2.8068 


0.2711 


9.498 


Alkanes 


2.014 






5.54 


0.025 





1.848 


1.026 






TABLE S3: Values of the constants used for the dielectric response e(aj) of Silicon H, Water and d. 



to analyze the fluctuations. In particular, infrared and microwave frequencies are very important as confirmed 
convincingly by Surface Force Apparatus experiments Q . 

An effective Hamaker "constant" is then defined by H {l2,T) = —12tt12Umi,M2 ih,T). For interlayer water thick- 
nesses dw = h + '^dhead larger than a nanometer, pairwise summation and Lifshitz theory give different values for the 
Hamaker constant (see Fig. IS7p . For our experimental values {d^ < 3 nm), the van der Waals interaction given by 
Lifshitz theory is well described by : -H{l2 = 0,T)/12tt {dyj + 2dheadf , with H{l2 = 0,T) ^ 5.3 x IQ-^i J. 



III. SELF-CONSISTENT THEORY (SC) AND EFFECTIVE POTENTIAL IN ASYMETRICAL 

POTENTIAL 

The renormalization of the microscopic interaction potential by the thermal fluctuations is a complex problem of 
modern statistical physic. Helfrich first realized that, in addition to the "direct" electrostatic, van der Waals and 
hydration forces @ , the long-range "effective" steric interaction generated by the thermal fiuctuations of confined fiex- 
ible membranes is an essential contribution to the total free energy of interaction Q • The effective potential, average 
bilayer position, and fiuctuation amplitude are thus coupled quantities which must be self-consistently determined. 
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o Lifshitz approach 
— Pairwise summation 

1271 U 



° 4 
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o Lifschitz approach 
— Pairwise summation 
12ji U 
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0.1 



I2 (nm) 




10 



FIG. S7: EfTective Hamaker constant vs I2 for a Si/Water/Bilayer/Water/Bilayer/ Water, pairwise summation (black solid line) 
and expression of the form —H{l2 = 0, r)/127rZ2 (black dashed line). In inset, zoom on the range of interest for our experiments. 

A. Effective potential theory 

Effective potential theories consist in adding an entropic (fluctuations) term C/fi (z) to the direct potential U {z). 
The dw equilibrium value is thus given by the minimum of U[^^- {dw) = {U + Ua)' {z = dyj) = 0. One can also compute 
the second derivative of the external potential at the equilibrium position C/j"^ {dw). 

• Pure hard wall interaction (hard confinement) was first considered in Q leading to an entropic contribution : 



(5) 



with ch ranging from 0.08 to 0.2 [1, Q in symmetrical case. 



• Podgornik and Parsegian have extended Helfrich's approach (l0| to take into account hydration repulsion and 
van der Waals attraction in the so-called "soft" potential [lH , 



Ur = Usoit = TTkBT/l&yJPh/K/zh cyi^{~ Z / nZh) ■ 



(6) 



B. Self-consistent theory 

Self-consistent theory developed by considers the position fluctuations of a membrane close to a substrate in 
an external potential. It is an alternative approach to compute the average membrane position and its root mean 
square fluctuation amplitude in the range of moderate fluctuations much below the unbinding transition. For arbitrary 
potentials, partition function mean membrane-substrate distance d^, fluctuation amplitude r.m.s. uth and mean 
value of second derivative of the external potential can be expressed as self-consistent equations : 



13 



Z = 



'th 



iu") 



( 16Kcrj\ 3 (z - d^)' 

dzcxp -^U {z) - 



[kBTY 



(z) = — / dzz cxp ^JJ (z) „ „ — — 



't/i. 



/ dziz-d^fe^p(-^^Uiz) 
^ J \ {kBT) 

1 /" , rr//^ N / ISkct^ 3(2;-d.u,)^\ 
— / dzC/ (z) cxp ^C/ (z) — — ■ 



'th 



3 (z - dw) 



C. Comparison between both approaches 

We have solved the self-consistent equations with an asymmetric external potential U (z) composed of hydration 
pressure term I7hyd(2^) = J^/t^/i exp(— z/z/j) and an osmotic pressure C/osmo(-z) = Pz (see Fig. [S8] inset). By varying 
P/, and z;i from (P/j = 0.5 • 10^ Pa, z/, = 0.18 nm) to [P^ = 0.5 • 10^ Pa, z/j = 0.005 nm) wc were able to mimic a 
transition from a soft hydration repulsion to hard wall stcric repulsion. 



E 




1 2 

dw (nm) 



FIG. S8: Mean values of the external potential second derivative calculated using self-consistent approach (points) and using 
effective potential theory (lines) for soft-repulsion (□ and dashed line {Pn = 0.5 • 10^ Pa, zn = 0.18 nm)) and for hard repulsion 
(• and solid line (Pn = 0.5 • lO" Pa, Zh ~ 0.005 nm)). The inset shows the associated external potential. 



By comparing the results obtained by the self-consistent and "effective potential" theories fFig |S8[) wc observe that : 

• In the case of a soft repulsion, self-consistent theory (□) and effective potential theory (blue solid line) are in 
good agreement, in particular in the small fluctuation regime, and for various bending modulus ranging from 5 to 
200kBT (not shown here). It is impossible to adjust a "hard wall" model , whatever the bending modulus used. 
However, a perfect fit can be obtained with a "soft potential" ?7hyd + C^oft, where the only fit parameter is the 
screening coefficient n in the confinement potential t4oft(-z) oc exp(— z/riz/j). The value obtained for ri is n = 2.4, 
in very good agreement with the literature [l2j . Nevertheless, the self-consistent approach is underestimating 
the mean distance d^ and overestimating {U"). 
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• In the case of hard wall repulsion (•), it is only possible to fit the results with Helfrich effective potential: 
f^Hei = ch I niksT / zY , with ch ^ 0.5. Whatever the parameters used, a "soft potential" model fails to describe 
this hard wall system. The value of ch is twice as large as the values obtained in the symmetrical case. 

These results show that self-consistent theory is able to describe the continuous transition from a soft hydration 
repulsion to an hard steric one. 
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